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I.  INTRODUCTION 


There  are  a  number  of  relatively  new  satellite-based  techniques  used  to 
measure  various  components  of  the  thermospheric  neutral  wind  field.  A  Fabry- 
Perot  interferometer  (FPI)  (Ref.  1)  was  flown  on  Dynamics  Explorer-B  (DE-B) 
(Ref.  2)  for  remote  sensing  of  the  meridional  wind  component.  This  instrument 
measured  the  Doppler  shift  of  the  thermospheric  0(^D)  line  to  derive  winds  for 
altitudes  between  150  and  350  km  (Ref.  3)-  Subject  to  the  availability  of 
appropriate  emission  features,  this  technique  is  accurate  (30  m  s-1)  and  is 
presumably  adaptable  to  cross-track  wind  measurements.  The  satellite  electro¬ 
static  triaxial  accelerometer  (SETA)  has  also  been  used  to  measure  cross-track 
(horizontal)  winds  in  situ  (Ref.  4).  This  technique  is  applicable  in  the 
altitude  region  of  significant  satellite  drag,  usually  below  230  km,  and  its 
accuracy  is  dependent  on  the  mass  and  cross  section  of  the  vehicle  and  on  such 
factors  as  satellite  noise  and  vibration.  A  random  error  of  60  m  s-"1  at  170 
km  was  quoted  in  the  above  work. 

Techniques  for  wind  measurement  employing  neutral -mass  spectrometers  have 
also  been  demonstrated.  A  system  which  used  a  single  moving  baffle  to  modu¬ 
late  the  incoming  gas  stream  was  flown  by  Spencer  et  al .  (Ref.  5)  on  Atmo¬ 
sphere  Explorers  (AE)  -C,  -D,  and  -E.  The  vertical  cross-track  wind  component 
was  measured  on  AE-C  (Ref.  6)  while  the  horizontal  cross-track  (meridional) 
wind  component  was  measured  on  AE-E  with  50  m  s-1  or  better  accuracy  at  260  km 
altitude  (Ref.  7).  Exploiting  spin  modulation  of  the  fly-through  mode  signal, 
Knutson  et  al .  (Ref.  8)  used  the  Open-Source  Neutral-Mass  Spectrometer  (Refs. 
9,  10)  on  AE-D  to  measure  the  vertical  wind  component  at  auroral  latitudes. 
In  order  to  measure  both  cross-track  components  of  the  neutral  wind,  a  dual 
baffle  system  was  employed  by  the  wind  and  temperature  spectrometer  (WATS) 
(Ref.  11)  which  was  flown  on  DE-B. 

To  date,  there  are  no  published  reports  of  a  successful  in  situ  measure¬ 
ment  of  the  in-track  component  of  neutral  wind.  In  principle,  this  wind 
component  could  be  obtained  from  an  energy  analysis  of  the  incident  gas 
stream.  Philbrick  et  al.  (Ref.  12)  obtained  a  retarding  potential  curve  for 
neutrals  with  a  velocity  mass  spectrometer.  Hanson  and  Heelis  (Ref.  13) 
described  a  variety  of  techniques  for  measurement  of  ion  and  neutral  winds, 
including  the  one  for  neutral  retarding  potential  analysis  (RPA).  Since  then, 
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there  have  been  attempts  to  make  this  measurement  on  AE-E  (Ref.  13),  Pioneer 
Venus  Orbiter  (Ref.  14),  and  DE-B  (Ref.  11).  Sensing  of  the  in-track  wind 
component  is  required  to  measure  meridional  wind  at  all  latitudes  from  a 
polar-orbiting  satellite.  Combined  with  the  WATS  dual  baffle  technique,  this 
additional  capability  would  enable  a  complete  specif ication  of  the  neutral 
wind  as  was  done  on  DE-B  with  WATS  and  FPI  (Ref.  3).  However,  measurements  of 
all  wind  components  by  mass  spectrometer  would  be  obtained  entirely  in  situ 
and  would  be  independent  of  atmospheric  emission  features. 

The  upper  atmosphere  composition  spectrometer  (UACS)  was  flown  on  Depart¬ 
ment  of  Defense  Space  Test  Program  mission  S85-1  in  the  summer  of  1984.  The 
orbit  was  nearly  polar  at  ~  84  degree  inclination  with  equatorial  crossings  at 
approximately  10.5  and  22.5  hours  local  solar  time.  UACS  measurements  were 
obtained  at  various  times  from  July  19  through  August  7,  while  perigee  and 
apogee  were  202  and  253  km,  respectively. 

In  an  earlier  paper,  Kayser  et  al .  (Ref.  15)  described  the  UACS  in  some 
detail  and  included  preliminary  examples  of  normal  and  fly-through  mass 
spectra.  Also  shown  was  an  RPA  sweep  on  the  molecular  nitrogen  peak.  Since 
then,  an  analysis  of  the  RPA  mode  and  its  operation  has  been  completed.  The 
wind  data  obtained  from  this  mode  have  been  reduced  and  compared  with  previous 
wind  measurements  from  other  sources.  Although  the  UACS  RPA  data  set  is 
small,  it  is  sufficient  to  demonstrate  successful  application  of  the  technique 
to  neutral  wind  measurements. 


II.  INSTRUMENTATION 


Details  on  the  operation  of  UACS  were  covered  in  an  earlier  paper  (Ref. 
15).  A  brief  description  is  included  herein.  Figure  1  is  a  cross-sectional 
view  of  the  sensor  with  its  vacuum  cover  removed.  A  portion  of  the  cover 
mechanism  (A)  remains  but  is  out  of  view  of  the  sensor  due  to  collimator  B. 
The  incident  gas  enters  the  collimator  from  the  left  and  impinges  on  the  ion 
source  chamber  (C).  Two  redundant  electron  guns  (D),  which  are  normal  to  the 
sensor  axis  (and  to  each  other),  produce  a  50-vA  electron  beam  which  is  used 
to  ionize  the  gas.  The  ions  are  extracted  (normal  mode)  or  drift  (in  RPA  and 
fly-through  modes)  into  the  lens  (E)  which  accelerates  and  focuses  the  beam 
into  the  electric  quadrupole  filter  (F).  Ions  of  the  selected  mass  exit  the 
quadrupole  where  96(&  are  collected  in  the  electrometer  cup  (G).  The  remaining 
4J  exit  through  the  back  of  the  cup  and  are  deflected  by  the  electrostatic 
deflector  H  into  the  spiraltron  electron  multiplier  (I).  At  the  beginning  of 
Aeoh  0.1  second  sample  period,  a  decision  is  made  by  the  electrometer  monitor 
circuit  whether  to  activate  the  deflector.  In  this  way,  a  large  dynamic  range 
can  be  achieved  while  protecting  the  multiplier  from  excessive  count  rates. 
The  RPA  signal  is  proportional  to  the  ambient,  rather  than  source,  gas  density 
and  therefore  is  inherently  small.  At  the  given  high  altitude  of  operation, 
all  RPA  signals  were  measured  on  the  multiplier. 

Figure  2  shows  the  ionizing  chamber  in  greater  detail.  The  parts  shown 
are  all  cyl indr ical ly  symmetric  about  the  neutral  beam  direction.  As  in  Fig. 
1,  neutrals  enter  from  the  left,  are  ionized  by  the  crossed  electron  beam,  and 
are  extracted  into  the  lens  on  the  right.  The  75  eV  electron  beam  is  col¬ 
lected  on  the  circular  collector  ring  which  is  part  of  the  filament  control 
circuit.  For  RPA  work,  shells  1  and  2  are  driven  equally  at  a  voltage  which 
brings  the  average  exit  energy  (and  quadrupole  drift  energy)  of  the  selected 
species  to  10  eV.  Hence,  for  N£  which  impinges  with  8.9  eV  energy,  this 
voltage  is  1.1  V.  For  atomic  oxygen  the  voltage  increment  is  4.7  V.  With  the 
mass  filter  tuned  to  a  selected  mass,  the  retarding  grid  is  driven  in  steps  of 
0.25  V  between  =  6  and  16  V,  bracketing  the  expected  energy  distribution. 


The  sensor  output  signal  is  proportional  to  the  number  of  ions  with 
energies  exceeding  the  barrier  energy.  A  typical  signal  from  the  RPA  analysis 
at  233  kni  is  shown  in  Fig.  3-  Forty-one  voltage  steps  are  executed  at  a  rate 
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A  Schematic  of  the  UACS  Sensor.  The  ambient  gas  stream  enters  from 
the  left  through  collimator  B.  Ion  source  region  is  C,  D,  and  E.  1 
is  the  electric  quadrupole  mass  filter.  Components  of  the  dual 
detector  are  found  at  G,  H,  and  I. 
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Detail  of  the  Ion  Source  Ionization  Region.  Neutrals  enter  from  the 
left,  are  ionized  by  the  crossed  electron  beam  and  extracted  to  the 
right.  Shells  1  and  2  define  a  drift  region  in  the  RPA  mode  and  an 
extraction  region  in  normal  mode.  The  retarding  grid  is  held  at 
ground  for  normal  mode  but  is  driven  by  a  ramp  generator  to  perform 
the  RPA. 
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of  ten  steps  per  second.  Not  shown  are  normal-mode  point?  taken  before  and 
after  the  sweep  and  between  voltage  steps  20  and  21.  Also  not  shown  are 
points  in  the  tail  above  14  V,  that  are  not  significant  for  the  analysis,  and 
a  number  of  "housekeeping"  steps  that  are  executed  at  the  beginning  of  the 
timing  cycle.  With  extraneous  steps  included,  a  total  of  5  seconds  (39  km 
distance  along  the  track)  are  used  in  one  energy  measurement. 


III.  TECHNIQUE  AND  DATA  ANALYSIS 


In  an  ideal  system,  the  RPA  response  curve  would  be  proportional  to  an 
integral  of  che  incident  gas  energy  distribution  function,  f(E),  where  the 
integral  is  taken  over  all  energies  E=eV  (V>Vr)  and  where  Vr  is  the  instanta¬ 
neous  retarding  voltage.  Absolute  particle  energies  would  be  obtained  direct¬ 
ly  from  the  calibrated  RPA  sweep  voltage.  In  practice,  a  number  of  instru¬ 
mental  characteristics  modify  the  overall  transmission  factor.  These  charac¬ 
teristics  include  the  response  of  the  mass  filter,  any  field  leakage  from  the 
ion  lens  through  the  retarding  grid,  off-axis  ion  discrimination,  residual 
electric  fields  in  the  source  drift  region,  and  any  electron  beam  potential. 
In  consequence,  the  RPA  curve  is  proportional  to  an  integral  over  q(V,Vr,...) 
f(eV),  where  q  is  a  function  of  instrument-related  variables.  Due  to  the 
factors  represented  by  q,  the  inherent  absolute  calibration  is  lost  and  some 
laboratory  calibration  is  necessary. 

A  typical  UACS  RPA  curve  is  shown  in  Fig.  3.  Provided  that  q  is  a  slowly 
varying  function  of  V  over  the  energy  range  of  6  to  16  eV ,  the  main  effect  of 
small  wind-related  shifts  in  the  energy  of  the  incident  gas  is  to  change  the 
position  (and  to  a  lesser  degree,  the  shape  and  amplitude)  of  the  RPA  curve, 
in  this  circumstance,  changes  in  position  should  be  a  direct  measure  of  wind 
over  the  expected  (<1000  m/s)  range.  The  kinetic  energy  represented  by  the 
50 1  amplitude  point  is  eVp-Eg-dE  which  is,  to  a  first  approximation,  equiva¬ 
lent  to  the  velocity  expression  1/2  m(vs+w)2,  where  Vp  is  the  retarding  volt¬ 
age  at  the  50%  amplitude  point  on  the  RPA  curve,  Eq  is  any  energy  increment 
imparted  to  the  ion  by  the  source  field  configuration  (e-beam  or  drift  space 
related),  dE  is  a  gas  pressure  correction  to  Eq,  m  is  the  ion  mass,  and  vg  is 
the  satellite  velocity.  Solving  for  the  in-track  wind  w  gives 

w  .  [2(eVR-E0-dE)/m]1/2-v3.  ( 1) 

As  a  final  stec,  w  must  be  corrected  for  any  component  of  the  rotating  atmos¬ 
phere  parallel  to  the  satellite  velocity  vector. 

Preferably,  the  RPA  mode  would  be  calibrated  in  a  neutral  beam  facility 
to  establish  an  exact  relationship  between  the  RPA  voltage  scale  and  the 
neutral  beam  velocity.  Lacking  this,  a  rough  laboratory  calibration  was 


established  as  follows:  The  instrument  was  operated  with  a  static  sample  gas 
and  with  source  voltage  settings  which  would  produce  nominal  10  eV  ions  at  the 
retarding  grid.  Then  the  normal  mode  voltage  differential  between  shells  1 
and  2  was  reduced  in  steps  and  at  each  step  the  RPA  curve  was  measured.  From 
these  data,  energy  shifts  in  the  curve  due  to  combined  electron  beam  and 
source  potentials  could  be  calculated.  The  energy  shift  inferred  for  zero 
potential  difference  between  shells  was  taken  as  Eq.  To  this  laboratory 
value,  a  correction  was  added  to  account  for  adjustments  of  source  voltages 
made  during  preflight  calibration.  Eq  was  determined  to  be  +0.9  (+  0.2)  eV. 
Thus,  due  to  the  particular  combination  of  source  and  beam  potentials  in  the 
present  design,  the  position  of  the  RPA  curve  was  shifted  by  0.9  eV  toward 
higher  energies.  The  net  increment  is  positive  in  spite  of  the  negative  well 
effect  of  the  electron  beam  because  the  space  between  shells  1  and  2  allows 
the  positive  collector  ring  voltage  to  protrude  into  the  drift  space  along  the 
axis.  For  the  same  reason,  when  shells  1  and  2  are  driven  to  equal  voltages, 
there  remains  a  nonuniformity  in  the  field  of  the  drift  space,  which  appears 
as  an  energy  dispersion  in  the  exiting  ion  beam. 

In  flight,  for  the  purpose  of  measuring  shifts  in  the  retarding  curve 
(hence  wind),  it  would  suffice  to  fit  some  empirical  function  to  the  sweep 
data.  It  is  more  instructive,  however,  to  approximate  the  shape  of  the  curve 
from  the  ambient  gas  distribution  function  and  estimates  of  the  instrumental 
response  function  q.  The  laboratory-measured  mass-filter  response  to  changes 
in  ion  energy  was  approximated  by  a  capacitor  charging  curve  whose  amplitude 
rises  to  63%  at  A  eV  and  to  essentially  100%  at  10  eV.  Placement  of  dual 
high-transmission  grids  shown  in  Fig.  2  effectively  eliminates  leakage  of  the 
ion  lens  field  through  the  retarding  grid  structure.  Thus,  transmission  of 
ions  was  set  equal  to  unity  for  particle  energies  exceeding  eVp,  and  to  zero 
otherwise.  Off-axis  ion  transmission  was  set  equal  to  unity  for  any  ion 
passing  through  the  aperture  of  ion  lens  element,  and  to  zero  otherwise.  It 
is  likely  that  the  effective  radius  is  actually  smaller  than  the  aperture 
radius,  due  to  the  low  transmission  factor  for  ions  entering  the  mass  filter 
at  large  angles.  However,  using  this  radius  as  the  approximate  extent  of  the 
effective  ion-producing  volume  in  the  electron  beam,  and  employing  the 
calculated  source  electrostatic  field,  an  energy  dispersion  can  be  estimated 
for  the  ions.  As  derived  above,  the  average  potential  in  the  production 
volume  is  the  shell  voltage  plus  E  . 
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With  these  estimates  for  q,  a  trial  RPA  curve  can  be  calculated  for  the 
incident  gas  energy  distribution  function.  A  family  of  such  curves  represent¬ 
ing  different  4^  is  generated  and  compared  to  the  RPA  sweep  data.  The  curve 
with  the  best  fit  by  least-squares  criteria  yields  the  best  estimate  of  Vp. 
Such  a  result  is  illustrated  in  Fig.  3,  where  the  best-fit  curve  and  curves 
displaced  by  +  0.1  V  are  shown.  Even  though  some  of  the  approximations  used 
in  the  generation  of  the  curve  are  crude,  the  resulting  shape  is  a  satisfacto¬ 
ry  reproduction  of  the  data.  With  Vp  known,  Eq.  1  is  employed  to  find  w.  The 
least-squares  process  also  yields  an  error  estimate  for  w.  For  the  wind 
measurements  to  be  presented,  these  errors  are  typically  0.036  eV  (15.3  m  s” ^ ) 
at  201  km  altitude  and  0.049  eV  (20.9  m  s” 1 )  at  226  km. 

While  in  principle  any  of  the  peaks,  N+ ,  0+ ,  or  ®2*  cou^d  be  usec* 
for  the  RPA  technique,  all  of  the  S85-1  wind  measurements  presented  here  were 
derived  from  retarding  potential  analysis  of  the  major  singly-ionized  molecu¬ 
lar  nitrogen  peak  at  mass  28,  thus  avoiding  possible  complications  due  to 
contamination  by  dissociative  ionization.  Potential  contamination  of  the 
signal  by  ambient  ions  is  minimized  by  the  overall  ion -optical  design  of  the 
sensor  and  by  the  electric  field  configuration  upstream  of  the  ionization 
chamber . 

In  Fig.  4,  the  derived  in-track  wind  for  all  local  times  has  been  plotted 
versus  latitude  using  the  satellite  reference  frame.  A  wind  of  positive 
magnitude  is  a  head  wind.  In  this  figure  only,  wind  was  calculated  with 
dE=0.  Turee  separate  data  segments  are  shown.  Data  represented  by  (-)  are 
all  part  of  a  southbound  10.5  LST  pass.  High  northern  latitude  data  (  +  )  are 
from  a  second  southbound  pass.  Points  (+)  between  -30  and  +12  degrees  are 
from  a  north-bound  pass  at  22.5  LST.  Due  to  the  higher  altitude  of  the  north¬ 
bound  pass,  the  signal  was  weaker  and  the  data  scatter  larger  than  for  the 
lower  altitude  measurements. 

A  statistical  analysis  of  the  data  in  Fig.  4  indicates  that  the  errors 
deduced  from  the  above  least-squares  process  should  be  multiplied  by  a  factor 
of  approximately  1.32  to  achieve  consistency  with  overall  along-the-track 
scatter.  If  this  factor-  is  attributed  entirely  to  technique,  then  typical 
measurement  errors  of  20  and  28  m  s'1  are  calculated  for  the  low-  and  high- 
altitude  passes  respectively.  Because  of  the  finite  angular  sensitivity  of 
the  sensor  and  the  sharply  peaked  entrance  angle  of  the  neutral  beam,  a 


NORMALIZED  SIGNAL 


Fig.  3.  Normalized  Detector  Signal  (Data)  for  N2  vs  Retarding  Potential 

Voltage.  A  best  fit  curve  is  shown  along  with  two  curves  which  are 
displaced  0.1  V  either  side  of  the  best  fit  voltage.  The  curves  are 
the  neutral  gas  distribution  function  convoluted  with  the  instrument 
response  function. 
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constant  cross  wind  will  change  the  amplitude  of  the  RPA  curve.  As  a  conse¬ 
quence,  variations  in  the  cross  wind  on  the  time  scale  of  a  second  could 
modulate  the  instantaneous  RPA  signal,  introducing  a  noise  component.  This 
seemingly  random  variation  would  contribute  to  the  measurement  error.  It  is 
likely,  therefore,  that  the  factor  of  1.32  represents  some  combination  of 
technique-related  errors  such  as  cross-track  wind  modulation  and  actual  short¬ 
term  variation  in  the  in-track  wind  component. 

To  the  extent  that  laboratory  data  can  be  used  to  fix  the  quantities  in 
Eq.  1,  the  measurement  technique  has  no  free  parameters.  any  change  in  the 
absolute  calibration  of  the  RPA  energy  scale  will  be  seen  in  Fig.  4  as  a 
uniform  vertical  displacement  of  all  points.  A  0.1  V  change  in  calibration 
would  shift  the  pattern  by  43  m  s”\  It  will  be  demonstrated  in  succeeding 
sections  that  changes  much  larger  than  this  would  lead  to  systematic  disagree¬ 
ments  with  existing  data.  Calibration  is  also  influenced  by  one  other  parame¬ 
ter,  namely  the  pressure-related  energy  increment,  dE,  induced  by  the  electron 
beam. 

Linder  and  Hernquist  (Ref.  16)  give  an  expression  for  the  relationship 
between  neutral  gas  pressure  and  the  depth  of  the  potential  well  in  an  elec¬ 
tron  beam.  Their  expression  can  be  used  to  derive  a  simplified  estimate  for 
changes  in  potential  with  pressure,  giving 

dE  =  Ew  ln(p2/p, )  (2) 

where  Ew  is  the  kinetic  energy  of  the  escaping  ions  at  the  walls  of  the 
chamber  and  p2  and  p1  are  the  final  and  initial  pressures,  respectively.  A 
value  of  0.112  eV  was  chosen  for  Ew.  This  is  four  times  the  energy  of  thermal 
ions  at  300  K  but  is  consistent  with  the  potential  well  depth  inferred  from 
lab  data.  p^  corresponds  to  the  pressure  in  laboratory  measurements  wherein 
the  value  for  Eq  was  inferred.  It  also  corresponds  to  low  altitude  near  the 
30-degree  latitude  point  in  Fig.  4.  Thus,  the  greatest  changes  from  the  dE 
correction  are  in  the  higher  altitude  22.5  LST  data,  shifting  this  segment 
by  -  60  m  s'1  toward  zero.  In  Figs.  5  and  6,  these  corrected  data  have  been 
plotted  (-)  after  averaging  in  two-degree  bins.  In  this  form,  the  data  repre¬ 
sent  the  meridional  wind  component  except  at  latitudes  above  55  degrees,  where 
the  satellite  path  introduces  a  small  and  unknown  zonal  component  in  the  in- 
track  wind  velocity. 
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ALONG  THE  TRACK  WIND 
IN  SATELLITE  REFERENCE  FRAME 


LATITUDE  (deg) 


UACS  Uncorrected  Wind  Measurements  in  the  Satellite  Reference  Frame. 
Winds  at  high  northern  latitudes  and  winds  with  positive  values  near 
the  equator  were  measured  on  southbound  passes  at  10.5  LST  and 
205(+5)  km.  Winds  at  low  latitudes  plotted  with  +  symbols  were 
measured  on  a  northbound  pass  at  22.5  LST  and  228(+7)  km.  In  the 
satellite  frame  positive  wind  is  a  headwind. 
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Combined  Data  Sets  of  the  Summer  Nighttime  Meridional  Wind.  Dashes 
are  the  corrected  and  smoothed  UACS  measurements.  Letters  and 
letters  with  vertical  bars  are  previously  published  data  from  AE-E 
and  a  variety  of  ground  stations.  As  much  as  possible,  the  data 
shown  are  associated  with  solar  minimum  condition.  Winds  to  the 
north  are  positive. 
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LST  - 10.5  hr 


MERIDIONAL 
NEUTRAL  WIND 


LATITUDE  (deg) 


Combined  Data  Sets  of  the  Sultrier  Daytime  Meridional  Wind.  Dashes  are 
the  corrected  and  smoothed  UACS  measurements.  Letters  and  letters 
with  vertical  bars  are  previously  published  data  from  AE-E  and  a 
variety  of  ground  stations.  As  much  as  possible,  the  data  shown  are 
associated  with  solar  minimum  condition.  Winds  to  the  north  are 
positive. 


The  measured  data  should  be  most  typical  of  solar  minimum  and  geomagneti- 
cally  quiet  times.  Figure  7  shows  both  the  solar  flux  at  1  AU  (Sa)  and  daily 


magnetic  index  (Ap)  for  the  latter  part  of  July  1984.  On  the  two  days  when 


winds  were  measured  (vertical  arrows),  the  solar  flux  remained  near  85  and  the 
magnetic  index  was  below  20.  This  circumstance  facilitated  a  comparison  with 
quiet-time  data  obtained  from  other  sources.  In  the  next  section  details  of 
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IV.  A  COMPARISON  WITH  OTHER  MEASUREMENTS  OF  THE  MERIDIONAL  WIND  COMPONENT 


The  present  results  can  be  compared  with  published  meridional  wind  meas¬ 
urements  which  were  derived  from  incoherent-scatter  radar  and  Fabry-Perot 
interferometer  techniques.  Fabry-Perot  measurements  are  available  only  for 
nighttime  hours,  but  at  a  variety  of  latitudes,  whereas  radar  data  are  avail¬ 
able  for  both  day  and  night,  but  with  limited  latitude  coverage.  Because  UACS 
data  were  obtained  in  a  period  of  low  solar  and  magnetic  activity,  the  compar¬ 
ison  excludes  much  recent  ground-based  data  taken  during  the  present  solar 
cycle.  Table  1  summaries  the  selected  ground-based  data  set  for  local  times 
of  10.5  and  22.5  hours  that  should  most  closely  approximate  conditions  during 
the  flight  of  S85-1.  These  data  are  plotted  along  with  the  S85-1  data  in 
Figs.  5  and  6. 


The  table  lists  for  each  data  value,  the  experiment  year  and  day,  local 


solar  time,  solar  flux  at  10.7  cm,  daily  magnetic  index  A  ,  measurement  alti¬ 


tude  in  km,  method,  a  one  letter  index,  and  source  reference.  Some  measure¬ 
ments  were  originally  reported  as  monthly  or  seasonal  averages.  In  other 
cases,  a  range  of  wind  values  was  inferred  by  this  author  to  be  typical  of 
conditions  reported  in  the  original  reference.  Thus,  the  day  column  includes 
designations  such  as  Wr,  Sr,  and  Eq  for  Sumner,  Winter,  and  Equinox,  as  well 
as  sequences  of  single-letter  abbreviations  for  months  (M-J-J  is  May-June- 
July,  for  instance).  Single  data  values  may  represent  single  measurements  or 
averages.  A  pair  of  numbers  in  parentheses  represents  a  range  of  veloci¬ 
ties.  Positive  values  indicate  a  northward  wind. 


The  literature  includes  numerous  examples  of  enhanced  meridional  wind 

near  midnight  [see,  for  example,  Amayenc  (Ref.  17)  and  Burnside  (Ref.  1 8 )  ] , 

with  much  of  the  increase  after  20-21  hours.  Care  was  taken  in  the  prepara¬ 
tion  of  Table  1  to  choose  nighttime  values  as  close  to  22.5  hours  as 

possible.  In  some  cases  it  was  necessary  to  make  interpolations  or  small 

extrapolations  from  later  hours.  In  two  cases  values  at  midnight  are  listed. 
As  much  as  possible,  measurements  at  low  solar  flux  and  low  solar  activity 
were  used.  Ideally,  wind  velocities  should  be  corrected  to  the  S85-1  satel¬ 
lite  altitude,  but  there  is  no  reliable  way  to  do  this.  In  the  example  of 
Burnside  et  al .  (Ref.  30),  even  monotonicity  with  altitude  may  not  be  a  good 
assumption . 
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Table  1. 

Ground- 

-Based 

Measurements  of  the 

Meridional  Wind 

Experiment 
Year  Day 

LST 

(hr) 

F10.7 

Ap 

Altitude 

(km) 

Method 

ID 

Wind 
(m/s ) 

Reference 

73-75 

Wr 

22.5 

low 

low 

250 

F-P 

a 

(-75, o: 

XRef.  19) 

73-77 

Sr 

22.5 

low 

all 

-250 

F-P 

b 

-40 

(Ref.  20) 

Sr 

-70 

Wr 

low 

-50 

86 

Jan 

14-17 

22.5 

75 

4.5 

300 

radar 

c 

-90 

(Ref.  18) 

77 

Aug 

21 

22.5 

90.8 

5 

-250 

F-P 

d 

0 

(Ref.  21) 

20 

22.5 

88.7 

8 

-250 

F-P 

e 

20 

(Ref.  22) 

76-77 

Sr 

22.5 

low 

low 

250 

F-P 

f 

(-80, 

-35) 

(Ref.  23) 

76 

Wr 

22.5 

low 

low 

250 

F-P 

f 

( - 140, 
-100) 

(Ref.  23) 

72-73 

M-J-J 

22.0 

med 

low 

250 

F-P 

g 

15 

(Ref.  24) 

A-S-0 

90 

M-J-J 

0.0 

100 

A-S-0 

100 

N-D 

22.0 

(50, 

150) 

76 

Jan  2 

0.0 

73.6 

4 

250 

F-P 

h 

60 

(Ref.  25) 

22.5 

-175 

75 

Nov 

22.5 

low 

low 

250 

F-P 

i 

(-50, 

(Ref.  26) 

0) 

Dec 

20 

85 

May 

22.5 

low 

low 

300 

radar 

j 

-125 

(Ref.  27) 

71-78 

Sr 

1 1.0 

all 

all 

all 

radar 

k 

-60 

(Ref.  28) 

Eq 

0 

Wr 

70 
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Table  1.  Ground-Based  Measurements  of  the  Meridional  Wind  (Continued) 


Experiment 
Year  Day 

LST 

(hr) 

F10.7 

Ap 

Altitude 
( km) 

Method 

ID 

Wind 
(m/s ) 

71 

July 

1 1.0 

300 

radar 

1 

(Ref.  17) 

6 

122.8 

9 

35 

7 

1 18.5 

3 

62 

8 

1 12.9 

8 

42 

9 

104.2 

6 

93 

74 

Aug  1  1 

.0 

200 

radar 

m 

(Ref.  29) 

The  data  of  Table  1  are  plotted  in  Figs.  5  (22.5  LST)  and  6  (10.5  LST'i 
using  the  single  letter  identifiers.  Capital  letters  represent  data  plotted 
at  the  measurement  location,  whereas  lowercase  letters  represent  data  plotted 
in  the  hemisphere  opposite  from  the  measurement.  Thus,  "b"  of  Fig.  5  was 
measured  in  Northern  hemisphere  winter  and  has  been  plotted  here  as  a  Southern 
hemisphere  winter  value  with  a  sign  change  in  wind  direction.  Table  2  lists 
ground  station  latitudes  and  longitudes  along  with  the  single  letter  identi¬ 
fiers.  In  the  example  of  "b",  the  source  was  Fritz  Feak.  Other  measurements 
at  this  location  in  summer  were  reported  for  low  sky  elevation  angles,  giving 
winds  at  35  and  45  degrees  latitude.  Whore  a  range  of  values  is  reported,  the 
letter  identifier  is  found  centered  in  a  vertical  bracket.  No  attempt  has 
been  made  to  assign  error  bars  to  the  referenced  measurements. 


Table  3  is  extracted  from  meridional  wind  data  published  by  Spencer  et 
al .  (Ref.  7)  and  Wharton  et  al .  (Ref.  3D.  These  wind  measurements  were  all 
made  in  1977  from  the  Atmosphere  Explorer-E  (AE-E)  satellite  before  the  onset 
of  the  present  solar  cycle.  The  orbit  of  AE-E  was  nearly  equatorial  with  an 
inclination  of  -19  degrees  and  altitude  between  255  and  280  km.  The  data  of 
Table  3  are  treated  in  a  manner  similar  to  the  data  of  Table  1  which,  f'or 
these  particular  data  samples,  gives  the  concentration  of  points  plotted  in 
Fig.  5  near  the  latitude  of  Arecibo. 


Referring  to  Fig.  5,  the  S85-1  nighttime  measurements  between  latitudes 
-37  and  +12  degrees  are  in  good  agreement  with  the  meridional  wind  profile 
implied  by  the  ground-based  data  (points  c,  d,  e,  and  h)  and  are  bounded  by 
the  AE-E  points  t  and  w.  The  comparison  of  S85-1  data  for  10.5  LST  with  otner 
measurements  (Fig.  6)  is  also  reasonable,  but  tentative,  considering  the  small 
number  of  other  points  available. 


As  was  mentioned  earlier,  a  zonal  wind  component  may  be  present  in  the 
UACS  data  at  latitudes  above  55  degrees  N.  latitude,  placing  those  values  at 
the  high  (winter)  end  of  the  range  measured  at  Chatanika.  On  the  other  hand, 
as  explained  by  Wickwar  et  al.  (Ref.  28),  the  lower  end  of  the  range  was 
determined  by  a  very  few  summer  observations.  We  cannot  be  certain  how  rep¬ 
resentative  these  points  might  be.  The  most  curious  wind  values  obtained  by 
'JACS  are  those  indicating  northward  flow  at  10.5  LST  between  the  equator  and  8 
degrees  N.  latitude.  While  this  flow  direction  does  not  agree  with  a  simple 
model  of  flow  away  from  the  solar  bulge,  there  is  a  supporting  data  point  (r) 
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from  Wharton  et  aJ .  (Ref.  31).  Also,  in  Fig.  6  the  range  of  values  measured 
at  Arecibo  (M)  would  appear  to  agree  with  reasonable  extrapolations  of  the 
S85- 1  data,  both  from  the  south  and  north.  Overall,  the  above  comparisons  of 
UACS  wind  measurements  with  other  available  data  provide  a  high  degree  of 
confidence  that  the  RPA  technique  was  working  properly. 

Beyond  the  confirmation  of  an  instrumental  technique,  the  meridional  wind 
profiles  of  Figs.  5  and  6  provide  an  interesting  glance  at  thermospheric 
dynamics  for  typical  solar  minimum  conditions  at  -  250  km  and  22.5  LST.  In 
Fig.  5  the  overall  profile  begins  at  southern  winter  mid-latitudes,  with 
northerly  flow  toward  the  equator.  This  flow  diminishes  and  reverses  in  the 
vicinity  of  10  degrees  N.  At  the  latitude  of  Arecibo,  the  flow  is  moderate  to 
strong  but  in  the  equatorward  direction,  as  has  been  regularly  observed  (Ref. 
32).  The  equatorward  flow  continues  at  higher  latitudes,  but  with  diminishing 
velocity  above  40  degrees  N.  A  few  data  points  for  25-30  degrees  N  would  be 
very  useful.  In  their  absence  it  is  still  reasonable  to  infer  that  in  summer, 
at  this  local  time  and  low  solar  flux,  the  equatorward  flow  velocity  maximizes 
around  25  degrees  N  and  is  typically  100  m  s_1 . 

The  summer  flow  picture  agrees  well  with  a  number  of  general  observa¬ 
tions.  Burnside  et  al .  (Ref.  32)  note  that  in  sunnier  the  equatorward  meridio¬ 
nal  wind  is  stronger  to  the  north.  Their  estimate  for  the  wind  gradient  of 
50-100  m  s'1  per  1000  km  agrees  well  with  a  value  of  60-80  m  s-1  inferred  from 
Fig.  5.  Also,  in  their  estimate,  the  meridional  wind  goes  to  zero  500  km  to 
the  south,  in  good  agreement  with  the  crossover  point  implied  by  the  figure. 
Some  caution  is  necessary  in  these  comparisons  because  the  wind  gradients  may 
vary  with  solar  activity  and  the  data  of  the  last  reference  were  obtained  at 
much  higher  values  of  F^  y . 

The  winter  hemisphere  wind  profile  may  be  more  controvers ial .  Burnside 
et  al .  (Ref.  30)  report  quiescent  or  even  poleward  flow  for  22.5  LST  at 
Arecibo  in  winter,  with  the  F10  y  solar  flux  near  200.  On  the  other  hand, 
point  C  in  Fig.  5  is  derived  from  Burnside  et  al.  (Ref.  18)  when  the  solar 
flux  was  75.  This  point  is  not  inconsistent  with  the  UACS  data  and  the 
Adelaide  data  (point  G).  The  winter  profile  therefore  may  be  typical  of  solar 
minimum  only  or  may  represent  one  extreme  of  wind  variability. 


A  final  point  of  interest  is  the  comparison  of  Arecibo  data  (Ref.  27) 
with  winds  from  AE-E  (Refs.  7,  31).  In  Fig.  5,  a  number  of  AE-E  points  fall 
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within  the  Arecibo  latitudes,  some  after  a  reflection  in  latitude  and 
season.  There  is  a  high  degree  of  consistency  among  these  data  and  reasonable 
agreement  with  the  ground-based  result  (J),  allowing  for  day-to-day 
variability.  Not  only  is  this  a  good  cross-calibration  between  techniques, 
but  it  offers  some  justification  for  the  use  of  hemispherically  "transplanted" 
data,  at  least  at  low  latitudes. 
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V.  DISCUSSION 


Overall  consistency  of  the  UACS  wind  measurements  with  previously  pub¬ 
lished  data  was  demonstrated  in  the  last  section.  In  a  sense,  it  was  fortu¬ 
nate  that  these  measurements  were  done  in  a  higher  than  anticipated  orbit 
because  the  bulk  of  comparison  data  had  been  obtained  at  these  altitudes.  On 
the  other  hand,  because  of  the  higher  altitude  of  S85  —  1 ,  there  are  several 
aspects  of  the  measurement  technique  that  were  not  fully  demonstrated  but 
which  are  important  for  lower  altitude  flights. 

Only  4%  of  the  available  count  rate  was  used  to  produce  the  RPA  signal. 
As  a  consequence,  the  random  error  of  a  single  4-second  wind  measurement  was 
on  the  order  of  15  to  21  m  s-^.  Had  the  instrument  been  designed  specifically 
for  wind  measurements  in  a  narrow  altitude  range,  a  factor  of  up  to  25  in¬ 
crease  in  signal  would  have  been  available.  While  the  implied  factor  of  5 
decrease  in  measurement  random  error  might  not  be  realized,  it  is  probably 
reasonable  to  anticipate  relative  errors  of  10-15  m  s” 1  for  a  repeat  experi¬ 
ment  which  is  better  matched  to  the  altitude  range.  Given  a  more  intense 
signal,  another  possibility  is  to  sample  at  a  higher  rate  to  achieve  in¬ 
creased  al ong-the-track  resolution. 

Thus,  w  i  *■  h  minor  changes  in  the  UACS  design,  it  is  expected  that  in-track 
wird  km  altitude  circular  could  be  measured  to  20  m  s-1  (relative 

error,  .er  16  km.  At  lower  altitudes,  signal  strength  improves  with  the 
scale  height  of  N2.  Absolute  calibration  of  the  instrument  for  wind  to  43  m 
s" 1  (0.1  Volt  for  N2)  or  better  now  apocars  feasible.  Given  a  large  database 
of  such  measurements,  it  may  also  be  possible  to  define  aposteriori,  the  zero 
of  the  RPA,  from  characteristics  of  the  observed  wind  field.  In  general,  the 
measurement  technique  presented  here  should  be  competitive  with  existing 
ground-  and  satellite-based  techniques. 

When  the  same  technique  is  used  on  lower-flying  satellites,  charge  neu¬ 
tralization  in  the  ionizing  electron  beam  may  produce  an  altitude-dependent 
offset  (dE)  in  the  measured  ion  energy  and  hence  in  the  inferred  wind  speed. 
This  effect  is  dependent  on  the  exact  geometry  and  configuration  of  the  ion 
chamber,  which  determine  the  degree  to  which  ions  produced  from  the  stagnated 
gas  component  are  trapped  in  the  region  of  the  electron  beam  (Ref.  33).  In  a 


/■/vV’ 


,  A  A  /.AAA  AAAA  A  AA  A  A  V 


,  •  .  «  .  A  A  ».  A  A  A  A  A  A  A  A  A  A  A  A  A  »  .  A  A  ^ 


design  with  efficient  trapping,  an  ion  created  from  an  incident  neutral  may 
see,  if  anything,  a  slight  net  positive  potential  at  the  electron  beam  and  be 
analyzed  as  having  a  "head-wind"  component  of  energy. 


At  the  other  extreme,  in  the  UACS  design  there  is  an  axial  electric  field 
component  in  the  ion  source  which  should  prevent  positive  ion  accumulation  at 
low  pressures.  In  this  case,  the  unshielded  electron  beam  imparts  a  net  nega¬ 
tive  energy  increment  to  incoming  particles.  In  combination  with  the  axial 
electric  field,  this  negative  well  broadens  the  energy  distribution  of  the 
incoming  gas  and  restricts  the  capability  to  measure  neutral  temperature. 

With  higher  pressures  at  lower  altitudes,  trapping  of  thermalized  ions 
should  resume  and  the  net  energy  increment  should  change  accordingly.  From 
laboratory  and  flight  data  it  appears  that,  without  corrections,  UACS  wind 
measurements  would  De  significantly  affected  at  an  altitude  of  150  km  with  an 
energy  offset  equivalent  to  ~  200  m  s_1  in  wind  speed.  Therefore,  a  source 
aesign  that  promotes  efficient  thermal  ion  trapping  over  the  operating  pres¬ 
sure  range  would  minimize  the  potential  well  effect  and  simplify  analysis  of 
wind  measurements.  The  shape  of  the  retarding  potential  curve  would  be  less 
altered,  leading  to  more  sensitive  temperature  measurements. 

Finally,  the  intercomparison  of  data,  done  here  to  validate  an  experi¬ 
mental  technique,  also  provides  an  example  of  the  potential  value  of  coordi¬ 
nated  ground-  and  satellite-based  measurements.  It  is  likely  that  the  com¬ 
bined  data  sets  in  Fig.  5  provide  a  better  view  of  the  low-  and  mid- 
latitudinal  meridional  wind  profile  near  midnight  than  has  been  previously 
published.  From  this  modest  exercise,  it  is  this  author's  impression  that 
existing  published  wind  data,  together  with  unpublished  data  that  must  exist 
from  ground  stations  and  satellites  AE-E  and  DE-B,  could  be  used  to  synthesize 
a  much  more  comprehensive  picture  of  average  wind  patterns  than  now  exists. 
This  can  certainly  be  done  with  new  measurements  by  satellite-borne  wind 
sensing  instruments,  such  as  UACS,  in  combination  with  the  ever-increasing 
number  of  IS-radars  and  F-P  observatories. 
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LABORATORY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  '‘architect-engineer”  for 
national  security  projects,  specializing  in  advanced  military  space  systems. 
Providing  research  support,  the  corporation’s  Laboratory  Operations  conducts 
experimental  and  theoretical  investigations  that  focus  on  the  application  of 
scientific  and  technical  advances  to  such  systems.  Vital  to  the  success  of 
these  investigations  is  the  technical  staff’s  wide-ranging  expertise  and  its 
ability  to  stay  current  with  new  developments.  This  expertise  is  enhanced  by 
a  research  program  aimed  at  dealing  with  the  many  problems  associated  with 
rapidly  evolving  space  systems.  Contributing  their  capabilities  to  the 
research  effort  are  these  individual  laboratories: 

Aerophyslcs  Laboratory:  Launch  vehicle  and  reentry  fluid  mechanics,  heat 
transfer  and  flight  dynamics;  chemical  and  electric  propulsion,  propellant 
chemistry,  chemical  dynamics,  environmental  chemistry,  trace  detection; 
spacecraft  structural  mechanics,  contamination,  thermal  and  structural 
control;  high  temperature  thermoraechanics ,  gas  kinetics  and  radiation;  cw  and 
pulsed  chemical  and  exciraer  laser  development  Including  chemical  kinetics, 
spectroscopy,  optical  resonators,  beam  control,  atmospheric  propagation,  laser 
effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions, 
atmospheric  optics,  light  scattering,  state-specific  chemical  reactions  and 
radiative  signatures  of  missile  plumes,  sensor  out-of -f ield-of -view  rejection, 
applied  laser  spectroscopy ,  laser  chemistry,  laser  optoelectronics,  solar  cell 
physics,  battery  electrochemistry,  space  vacuum  and  radiation  effects  on 
materials,  lubrication  and  surface  phenomena,  thermionic  emission,  photo¬ 
sensitive  materials  and  detectors,  atomic  frequency  standards,  and 
environmental  chemistry. 

Computer  Science  Laboratory;  Program  verification,  program  translation, 
performance-sensitive  system  design,  distributed  architectures  for  spaceborne 
computers,  fault-tolerant  computer  systems,  artificial  intelligence,  micro¬ 
electronics  applications,  communication  protocols,  and  computer  security. 

Electronics  Research  Laboratory:  Microelectronics,  solid-state  device 
physics,  compound  semiconductors,  radiation  hardening;  electro-optics ,  quantum 
electronics,  solid-state  lasers,  optical  propagation  and  communications; 
microwave  semiconductor  devices,  raicrowave/milliraeter  wave  measurements, 
diagnostics  and  radioraetry,  microwave/millimeter  wave  thermionic  devices; 
atomic  time  and  frequency  standards;  antennas,  rf  systems,  electromagnetic 
propagation  phenomena,  space  communication  systems. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metals, 
alloys,  ceramics,  polymers  and  their  composites,  and  new  forms  of  carbon;  non¬ 
destructive  evaluation,  component  failure  analysis  and  reliability;  fracture 
mechanics  and  stress  corrosion;  analysis  and  evaluation  of  materials  at 
cryogenic  and  elevated  temperatures  as  well  as  in  space  and  enemy-induced 
environments . 

Space  Sciences  Laboratory:  Magnetospherlc,  auroral  and  cosmic  ray 
physics,  wave-particle  interactions,  magnetospherlc  plasma  waves;  atmospheric 
and  ionospheric  physics,  density  and  composition  of  the  upper  atmosphere, 
remote  sensing  using  atmospheric  radiation;  solar  physics,  infrared  astronomy, 
infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth's  atmosphere,  ionosphere  and  magnetosphere; 
effects  of  electromagnetic  and  particulate  radiations  on  space  systems;  space 
instrumentation. 


i/vj 


vVC 


*  k*VI 


»•  J 


